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Abstract 

We propose a new reformulation of Yang-Mills theory in which three- and 
■ four-gluon self-interactions are eliminated at the price of introducing a suffi- 

ce , cient number of auxiliary fields. We discuss the validity of this reformulation 

\ in the possible applications such as dynamical gluon mass generation, color 

' confinement and glueball mass calculation. Moreover, we set up a new 1/iVc 

color expansion in the SU{Nc) Yang-Mills theory based on this reformulation. 
In fact, we give the Feynman rules of the 1/Nc expansion in the manifestly 
Lorentz covariant gauge. The Yang-Mills theory is defined on a non-trivial 
• vacuum where color-singlet transverse gluon pair condensations take place by 

D , the attractive gluonic self-interactions. This vacuum condensation provides a 

' common non-vanishing mass for all the gluons with color symmetry being pre- 

served. It is shown that the auxiliary fields become dynamical by acquiring the 
^ . kinetic term due to quantum corrections. Then the static potential between a 

H ' pair of color charges is derived as a combination of the Yukawa-type potential 

and the linear potential with non-vanishing string tension. The mass of the 
lightest scalar glueball is calculated as the ratio to the gluon mass. The ex- 
plicit calculations are performed as a partial resummation of the leading order 
diagrams for the small 't Hooft coupling. 
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1 Introduction 



The mechanism for providing the gauge boson with the mass in strong interactions 
is still under the debate, although fifty years are going to pass since the discovery of 
the non-Abelian gauge theory by Yang and Mills P . In this paper we argue that the 
normal massless vacuum of Yang-Mills theory is unstable, since a group singlet pair 
of transverse Yang-Mills gauge boson forms a bound state with a negative energy 
eigenvalue (Cooper pair) due to attractive force. Consequently, the stable vacuum 
is realized after the condensation of these pairs. The excitation spectrum has a gap 
without breaking the global gauge (color) symmetry for any coupling strength and 
the gauge bosons acquire a common mass. 

In the classical level, Yang-Mills theory is a scale-invariant (actually, conformal 
invariant) field theory describing massless gauge bosons. By the gauge principle, the 
existence of the mass term is prohibited in the classical Yang-Mills theory. In the 
quantum level, on the other hand, one has the possibility of the dynamical mass 
generation for the gauge boson as investigated long ago jHllZlEllSEllllEl • 

In the above context, among others, it is important to recall a pioneering work by 
Fukuda and Kugo j2] in which they discussed the transverse gluon pair condensation 
by using the variational technique for SU{N) Yang-Mills theory already in the late 
1970s. They adopted the Coulomb gauge in the Hamiltonian formalism and dealt 
with only the transverse modes with no ghost. The basic tools were the Cooper 
equation, Bogoliubov transformation and the gap equation. However, the resulting 
energy spectrum A(p) did not have the relativistically invariant form ^/J?~+ln?, since 
their analyses are based on the non-covariant formalism. In the second paper jS], 
Fukuda attacked the same problem in covariant form by making use of the effective 
action or potential and the resulting Schwinger- Dyson (SD) equation. He obtained 
qualitatively the same result as the first paper. These works demonstrated that the 
vacuum of the massive phase can have the lower energy than the normal vacuum. 




Figure 1: Gluon pair condensation in color singlet and spin singlet channel (center of 
mass frame). 

If two gauge particles form a bound state, the bound state should be a tachyon 
pole, since two gauge particles are massless. The existence of the tachyon bound 
state was also discussed by solving the Bethe-Salpeter (BS) equation in the ladder 
approximation jUEj where the Landau gauge A = is a preferred choice. From the 
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viewpoint of the dynamical symmetry breaking, the existence of Goldstone pole as a 
consequence of the condensation phenomena was discussed by several works [HIIZIIHI 
in the framework of SD and BS equations. 

However, it should be remarked that the solutions of the SD and BS equations 
strongly depend on the choice of gauge fixing parameter A. In fact, the result is 
strongly gauge dependent; for A > —1 the tachyon pole which is group singlet and 
spin singlet exists, whereas such a pole does not exist for A < —1. It is notorious that 
the systematic improvement is rather difficult in the SD and BS approach. 

In this paper, we discuss the problem based on a manifestly Lorentz covariant 
formulation. Therefore, we can obtain the manifestly relativistic result. Quite re- 
cently, it was shown ^2] that the transverse gluon pair condensation can be replaced 
by the on-shell BRST invariant vacuum condensate of mass dimension 2 which was 
proposed by the author ^QlE]- In this paper, we show that the l/N^ expansion can 
be set up to yield the condensation phenomena in the systematic way. The l/N^ 
expansion is a non-perturbative method which is systematically improvable, in con- 
trast to solving the approximate gap equation by adopting the ad hoc approximation. 
In our approach, the existence of the pair condensation is also a gauge independent 
phenomena in the sense; if it occurs at a certain gauge at all, it should occur in all 
the other gauges. 

The expansion was extensively used in the past works. The well-known exam- 
ples are among others the 0(N) symmetric four-fermion model (e.g., two-dimensional 
Gross-Neveu model) and the non-linear a model It is also remarkable that 

the three-dimensional four-fermion model become renormalizable in the expan- 
sion ^3], although it is non-renormalizable in weak-coupling perturbation expansion. 
The four-dimensional Yang-Mills theory in the l/N^ expansion might be related to 
the string theory jT3]. See e.g., JH] for reviews of large N expansion. 

It should be noted that the l/N^. expansion given in this paper is quite different 
from the conventional one for the Yang-Mills theory ^3^]. In this paper we adopt 
the Lorentz covariant gauge. The Maximal Abelian (MA) gauge will be treated in a 
separate paper [T7] . 

This paper is organized as follows. In section 2, we rewrite the original Yang-Mills 
action which includes the cubic and quartic gluon self-interactions to an equivalent 
theory without gluon self-interaction by introducing sufficient number of auxiliary 
fields. Consequently, the gluon field can be integrated out exactly. In section 3, we 
set up a new l/N^. expansion of the SU{N^ Yang-Mills theory on a non-trivial vacuum 
where the vacuum condensation of mass dimension two takes place. We give Feynman 
rules of the relevant 1/iVc expansion. The low-energy and high-energy behavior of the 
auxiliary field propagators are studied. As a by-product, the mass of the lightest 
glueball is given. In section 4, we obtain the effective potential of the theory and 
estimate the effective gluon mass by making use of the lattice data. In section 5, we 
calculate the static potential between a pair of color charges to the leading order of 
the 1/Nc expansion. We show that the static potential is composed of the Yukawa- 
type potential and the linear potential with the non- vanishing string tension which is 
proportional to the vacuum condensate of mass dimension 2. In section 6, we discuss 
the string representation of the Yang-Mills theory within the scheme of the l/N,, 
expansion. The final section is devoted to conclusion and discussion. Some of the 
technical materials are given in the Appendices A and B. 
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2 Reformulation based on the auxiliary field method 



2.1 Yang-Mills action 

The Euclidean version of the Yang-Mills Lagrangian is given by 

4 4 



+ X ■ X '^l^ (2-1) 

where J-'^^ is the non-Abelian field strength defined by JF^^, = df^A^ — d^A"^ + 
gjABCj^Bjc ^j^]^ A, 5, C = 1, ■ ■ ■ , dim(G) for the gauge group G. In what fol- 
lows, we use a notation := F ■ F, (F ■ G) := F^G^ = 2tr(FG), (F x G)^ := 
j!ABC jpBqc _ — 2ztr(T"^[F, G]) with the normalization for the generators T"^ of the 
Lie algebra W of the gauge group G: tr(T^T^) = ^S^^. Then the generating func- 
tional is given by the functional integral: 

Z[J] ■= J V^/^V^^V^^Vtf^ {~ / ^^^["^YM + -^^l} , (2.2) 

where we have included the source term = s^/^ ■ J^^ + ■ Jb + Jc ■ ^ + -h ■ ^ ■ 
Here the total Lagrangian of the Yang-Mills theory =^ym = -^ym + =^gf+fp, is given 
by adding the gauge fixing (GF) and Faddeev-Popov (FP) ghost term ^gf+fp to the 
pure Yang-Mills term =^ym where is the Nakanishi-Lautrup field, and and 
are the FP ghost and antighost fields. 

First, we introduce the second-rank antisymmetric tensor field B^y as the Hodge 
dual of the distorted field strength according to [THIIT^ I^ 

B^u ^ - <9,^^) + <yg{£/^ X O], (2.3) 

where p and a are two parameters to be determined later and the Hodge star (duality) 
operation * for the second-rank tensor is defined by *H^i, := \€fj_yp„H''" . For this 
purpose, we insert an identity into the functional integral ()2.2|) : 



Then the Lagrangian is modified as 

^YM =\b^, ■ B^" - ^B^" ■ *[p{dp^. - d,s^^) + ag{s^^ x s^,)] 



(2.4) 





4 ^ 


1 


- pa 




2 
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+ -^9{dpK - d,s^^) ■ (^^ X ^-) 

-/(^^ X O ■ (i^'' X .s^f""). (2.5) 
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For any values of the two parameters p and cr, the Lagrangian ()2.5|) is equivalent to 
the original Lagrangian ()2.1|) . Note that B^,^ is the auxiliary field without its kinetic 
term. 

Next, we rewrite the last term in ()2.5|1 which is quartic in the gluon field For 
SU{Nc), the quartic term has another expression^ 

K X o ■ i-^" X -^n ~ [K ■ -^ni-^u ■ -^n - k ■ ■ -^n] 

+ (i^ * ■ * - (i^ * O ■ (^/'' * i^^), (2.7) 

where we have defined a new product {Fi<G)^ := d^^^F^G^ = 2tr(T^{F, G}) with 
completely symmetric tensor d"^^*" for Nc > 3. For SU{2), we do not need to consider 
the last two terms in ()2.7|1 . since c?"^^"-^ = for SU{2), and hence the quartic term is 
greatly simplified 



X O ■ X j^'^) =(^4 ■ jy^)« ■ ^^) - (j4 ■ 0(^'^ ■ J^^). (2.8 
For later convenience, the identity ()2.7p is rewritten into 

X o . (^^ X ^-) =^(=< • =«/^)« ■ - ^s.^s^'^ 

+ ^(^^^^'^)2-T^,.T^-, (2.9 

by introducing the traceless symmetric tensors S^^ti, = S^^^ and = T^^ defined by 

1 
4 



S^u :=(■< ■ O - 7V(=^P ■ V = 0, (2.10) 



T^t : = (^^^0^-^(5^,(^,*^Y, TV = 0. (2.11) 

Moreover, we introduce the auxiliary scalar field and the auxiliary symmetric 
tensor fields G^^ for SU(2) and also ({>^ and for SU{Nc) {Nc > 3). For example, 
the auxiliary field (j) is introduced by inserting one more identity into ()2.2j) 

1 = y'D0exp|-|rf4x^(^0+^^^-^^)'|, (2.12a) 

where a^j, is an arbitrary parameter at this stage. Other auxiliary fields G^u,(^"^ and 
are introduced in the similar way. 



1 = j VG^^exp |- 1 d^x^ (g^, + , (2.12b) 

1 = |Py,^exp | - I d'x^ [ip^ + ^n^^'^ , (2.12c) 

1 = / ^^^,^.exp j- / d'x^ (V^i + It;^)'} • (2.12d) 



^We have used an identity for the structure constant (see [21] for a derivation): 

fABE fCDE ^ ^(S^CgBD _ gAD^BC^ ^ d-^^B^BD^ _ ^BCE^ADE^ (3.6) 

For the generator of the Lie algebra of the gauge group G, [T^,T^] = if^^'^T'^ and 
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Then the quartic interaction in ()2.5p is replaced as follows. 



2^2 



a^ + 3(l-(TV/^^c 



2 2 o 

+ '-^V,^ ■ V"" + ^V;„ ■ T-" + £il^^(i^^r,„ . T"". (2,13) 

Z Z o 

2.2 Lorentz gauge fixing and FP giiost 

We adopt the generalized Lorentz gauge or Curci-Ferrari (CF) gauge jl 
The gauge fixing (GF) and the Faddeev-Popov (FP) ghost term is given by 



^GF+FP = 5b5b ( ■ ■< + ■ ^) - f ^ ■ (2-14) 



= 5b (d^'i ■ + ■ - Kx'g"^ ■ X ^) j - ■ e^. (2.15) 

The conventional Lorentz covariant gauge corresponds to a' = 0. Performing the 
BRST transformation (Euclidean version), 

5B^^(a:) = ^^K]^(x) := d^^{x) + g{s^^{x) x ^(x)), (2.16a) 

5b^(x) = -l^(^(x) X ^(x)), (2.16b) 

5b<^(x) = -^(x), (2.16c) 

5g^(x) = 0, (2.16d) 
and anti-BRST transformation (Euclidean version), 

5b^^(x) = ^^K]^(x) := d^^{x) + ^(^^(x) X ^(x)), (2.17a) 

^B^(x) = -i(?(^(x) X ^»), (2.17b) 

^B^(x) = -^(x), (2.17c) 

5b^(x) = 0, (2.17d) 

where ^(x) = — ^(x) + g{^{x) x ^(x)), we obtain 

2" 



^GF+FP = - ■ ^ + ^ ■ + -^7(^ X ^") 



O- - - 

+ ^ ■ ^^^^K]'^ + -^'(^ X ^) ■ (^ X ^) (2.18) 



~2 



+ ^ ■ d^s^^ + — ^(^ X ^) 



+ ^-9^^^K]'^- ^(^'(^ X ^) ■ (^ X ^). (2.19) 
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The GF+FP term is further rewritten as 



5gf+fp= / d^x(^j^{d^^^)^ + {l-0'^-dM-^]'^ + ^'^-%Wf.'^ 



A 
2 



(2.20a) 



or 



5, 



GF+FP 



2A 



{d^s^,? + ■ d^d^^ + (1 - ■ (5m^ X ^) 



^(7=^/^ ■ X d^^) - -Ae(l - X ^) ■ (^ X '^) 

A 
2 



A-i9^^^-e^?(^x^))'| 
where we have defined the two parameters 

a'/2 ( 



A := a + a', ^ := 



a + a' 2A' 



(2.20b) 



(2.20c) 



In this form, it is easy to ehminate the Nakanishi-Lautrup field ^J^. 

For the Landau gauge A = 0, the ghost self-interaction disappears. In the gauges 
other than the Landau A 7^ 0, we treat the four ghost interaction in the similar way 
to the four gluon interaction: 



(^X ^) ■ (^X ^) = _(^.'^)2 + (^^^) . (^^^), 



(2.21; 



where we have used a fact that ^ ■ ^ = = ^ ■ and '^Tir^ = = '^*^ due 
to Grassmannian property. By introducing the auxiliary fields (j)FP and ippp through 
the identities: 



V(f)Fpexp 
Vcppp exp 



d X- 



2 

d'x^ 



pp + ■ 



ifpp + 



the four ghost interaction term is cast into the form: 



1 

~ 2 

(^FP ,2 



-A^(l - 0/(^ X ^) ■ (^ X ^) + 
a 



^FP ( , ^ , Opp 

pp -\- To ■ 10 



ypp 



2 

^^Ifpp ■lfFp + aFp(j)Fp'T^ ■ ^ + (TpplfiFP ■ -k'lo) 



(2.22a) 
(2.22b) 



aFp-\i{l-i)g^ — 



a'pp - Ae(l - 0/] (^"*^) ■ (^*'^). 

(2.23) 
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The auxiliary fields must be introduced such that the BRST symmetry is main- 
tained. For this purpose, the simplest choice of the BRST transformation for the 
auxiliary fields is given as 



= -6b{ ■ ) (2.24a) 



•1 

SbB^,, = -Sb Cipid^K - 9,^4) + agij^^ x <)]) , (2.24b) 

SbG^, = -Sb (^5^.) , (2.24c) 

5b^ = -Sb Q=< * , (2.24d) 

SbV^. = -Sb (^T^.) , (2.24e) 

(5B0i.p = -(5B(^-^), (2.24f) 



5b¥^fp = -5b ^*'^ . (2.24g) 



The original measure djj,{^) := T>^^{D^T>^T>l3S is known to be BRST invariant 
Therefore the additional measure djji{^) := VB^yV(j)VG^i,V(pW^iyV(f)FpV(pFP of the 
auxiliary fields is also BRST invariant, since it is easy to show that the Jacobian of 
the BRST transformation associated to the total measure (i/i($)(i/i(\E') is equal to 
one. In particular, we find 

(5b0(x) = -=<(x) ■ 9^^(x) = 9^i4(x) ■ ^(x) - d^'i^/^ix) ■ ^(x)]. (2.25) 

In the Landau gauge (9'^^(x) = 0, therefore, the BRST transformation of </'(x) is 
reduced to the total derivative and the spacetime average of (f){x) is BRST invariant, 
i.e., ^B / d^x(f){x) = 0. 

2.3 Total Yang-Mills action 

A purpose of this paper is to obtain an equivalent theory which does not involve the 
gluon fields ^ by integrating them out. The GF+FP term in the Lorentz gauge 
-^GF+FP (j2.2nb|l includes the linear and quadratic terms in the gluon field. 

The pure Yang-Mills Lagrangian is further rewritten as follows. For this purpose, 
we adjust the parameters to eliminate the cubic and quartic interaction terms in the 
field =2^. The cubic interaction term in (j2.5|l is eliminated by choosing 

p = a~\ (2.26) 

Moreover, by choosing the parameters as 

3{a' - 1)9' 4(1 -aV 
a^ = , (Tg = , = {3/2)[a - l)g , ay = 2[l - a )g , 

(2.27) 



the quartic gluon interaction terms in (|2.13p are also eliminated. Then a is only one 
undetermined parameter at this stage and it must be larger than one, i.e., o"^ > 1, to 
realize the condensation as shown later. 
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Due to the analysis 013123; the triple gluon interaction provides the attractive 
force, whereas the quartic gluon interaction can give the repulsive force among the 
gluons. Therefore, for the gluon pair condensation to occur, the attractive force 
coming from three-point vertex dominates the repulsive force from four-point vertex 
to yield the net attractive interaction between a pair of gluons. 

In our treatment, the original triple gluon vertex is completely eliminated by 
introducing the antisymmetric tensor B^i, and choosing p = cr~^. Therefore, the net 
attractive interaction is reflected in the relation cr^ > 1. In fact, if this is not the case, 
i.e., repulsive interaction 1 — o"^ > 0, the gluon pair condensation does not occur as 
will be shown later. 

Therefore, the pure Yang-Mills Lagrangian ()2.5p becomes at most quadratic 
in the gluon field and reduces to ^ 



+ ■ ^ + ■ (^^ ^ + ^V,, ■ V^'' + ^V,, ■ T^r (2.29) 



Moreover, the GF-I-FP term is reduced to 
1 

2A 



=^GF+FP =— + ■ dA"^ + (1 - OoK ■ {d^^ X ^) 



+ -y-^FP + -^^FP ■ <^FP + (TFP(f)Fp'^ ■ ^ + Cr'ppi^Fp ■ ^). 

(2.30) 

Here the four-ghost interaction terms in the GF-I-FP term are eliminated by choosing 

app = X^l - 0/, ^'fp = m - Og'- (2.31) 
Finally, we obtain the total Lagrangin which is at most quadratic in the gluon field 

I a2 , ^'fP 

+ -1^9fP + -l^^FP ■ ^FP 

+ <^ ■ (9^9^^ + aFp(t)Fp'i ■ <^ + a'pp^FP ■ (^* ^) 

+ ^^/^^^''^W,^ + ■ (2-32a) 



^If we want to leave two independent parameters p and cr, we have only to add 

+ ^—;^9{dt.< - 9.=0 • (^^ X ^'^). (2.28) 

Therefore, we can choose an interesting option p ~ a, i.e., B^^i, is the Hodge dual of -^^i,, although 
the cubic gluon interaction term remains in this case. 
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where we have defined J^i^^ and by 



:i-p'){d'6^,-d,d,)-x~%d. 



+ ctgS''^ (g,, - -6,.G/^ + ayrf^^^ (y;; - ^ V^p""") ' (2-32b) 
■■=J^ - ^pdu*B^, + (1 - 0^7(9,.^" X - X 9,^)^. (2.32c) 

Here is the source for the gluon field =2^. Note that is at most linear 

in the auxiliary field and contains the bilinear term in the ghost and antighost 
fields. The source term for the auxiliary field will be introduced later to avoid the 
complication. 

Thus, we can integrate out the gluon field exactly: 

+ —^<PfP + —^^FP ■ ^FP 

+ ^ ■ (9/,9^^ + CTfp^fP^ ■ ^ + a'pp^FP ■ (^* ^) 

- ^ A'^t'^^^'n'' A""} + ^lnDet[jr^^n- (2-33) 

When X = 0, (jp = = a'pp and the ghost self- interaction term disappears. Therefore, 
we do not need to introduce (ppp and ippp in ()2.33|) which decouple from the theory. 
Integrating out the gluon field and FP ghost fields we can obtain the 

effective theory written in terms of five kinds of auxiliary fields 5^^,, (p, (p"^, G^u, V^. 
The A 7^ case will be considered later. 
Some remarks are in order. 

1. The resulting theory has one undetermined parameter a. In the classical 
level, any value of a reproduces the original Yang-Mills theory, as is trivial from 
the above construction. However, in the quantum level, this is not the case. The 
renormalization urges the parameter a to run according to the renormalization scale, 
as will be discussed in Appendix A. 

2. For a special choice o"^ = 1, we have = cx^ = 1 and hence the coefficients 
vanish 0"$ = 0, which implies that all the auxiliary fields (j), G^i,, ip, V^y except for 
B^u decouple from the theory. Therefore, this case reproduces the field strength 
formulation [SI] of the Yang-Mills theory, which was also called the deformation of 
the topological BF theory [23122] with the Lagrangian: 

^BF = \B,y ■ B'^-' - '-B^- ■ *^^y. (2.34) 

(The signature of p and a can be absorbed by the redefinition of the fields ^ and 
Bf_iu-) When ^ 1, therefore, we have the scalar field 0, (f and the symmetric tensor 
fields G/jjy, Vfj_iy in addition to the antisymmetric tensor field S^j,. 

3. It is in principle possible to introduce the auxiliary field (pp corresponding to 
the square ^(-s^hys)^ of transverse gluon modes =2^hys (i-^-' ^i-^physi^) — 0) 

P^^exp |- / d'x^ (<j>P + ^«hys)')'| , (2.35) 
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where ^hys is both BRST invariant, (5B^hys(^) = O5 gauge invariant, d^^^p^y^^x) = 
0. Then the auxihary field (pT introduced in this way becomes BRST invariant, 
^b4>t{x) = 0, and gauge invariant, d^cpT^x) = 0. For example, such a transverse 
mode is written as 

<ys(^) = {S'' - d'^d^) = ^{.^ + g^l,^ + + ■ ■ ■ , (2.36) 

where can be determined order by order of the coupling constant g. The first 
two terms are $-|q-) = and = + ^[d^V(^i),V(^i)] where := ^dj^^, 

see section 5 and Appendix A of jHZI for more details. Needless to say, the transverse 
mode ()2.3(ij) reduces in the Abelian limit to the well-known expression: 

A^(x) = f 1 - V^V-) A{x) = A(x°, x) + ar / rf^y— ^ raM-'(x°, y). (2.37) 

V A / J 4:71 \x — y\ 

However, in this case, the gluon self-interaction can not be eliminated, since ^*hys is 
a non-local and furthermore non-linear quantity written in terms of in the non- 
Abelian gauge theory. (Even the all order expression of ()2.3fi|l is not yet known.) 
Therefore, the strategy adopted in this paper is difficult to be applied to this case 
which is more transparent from the physical point of view. 



2.4 Another reformulation 

We can introduce another antisymmetric tensor field B^^^' by 

- B^J := B^, - ip*{d,K - d^s^^)^, (2.38) 
which is equivalent to the identification as an auxiliary field: 

B^J ^ crg*{£/^ X ^,). (2.39) 
Then the Lagrangian is modified as 

^YM =\b^J ■ B^^' - '-B^^' ■ *[ag{^^ x <)] 
+ \9{d,K - d,^,) ■ {s^^ X ^'^) 

+ ^^^/(i^ X ■ {s^^ X s^""). (2.40) 

By this redefinition of the antisymmetric tensor field, the original cross term —^B^'^ ■ 
p*{d^s^v — dys/fj) has disappeared. The price is that the three-gluon interactions re- 
main in this case. The four- gluon interactions can be eliminated by applying the same 
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procedure as that given above. Thus we obtain the modified Yang-Mills Lagrangian^ 

In order to treat the p = a case, anyway, the three-gluon interaction can not be elimi- 
nated. The absence of the four-gluon interactions eliminates the two-loop diagrams in 
the Schwinger-Dyson equation for the full gluon propagator. Although we must deal 
simultaneously with a number of auxiliary field propagators and the cross propagator 
between the gluon and the antisymmetric tensor field, all the diagrams are of the one- 
loop type and there is no ambiguity encountered in treating the two-loop diagrams. 
If we set cr = 0, furthermore, the antisymmetric tensor field B^i^' decouples from the 
theory. Thus this reformulation helps to simplify the Schwinger-Dyson equation for 
the full gluon propagator, since the absence of cross propagator reduces the number 
of propagators to be investigated. The details will be given elsewhere pHj . 



3 Setting up the 1/Nc expansion 

In the previous section, we have rewritten the Yang-Mills theory into an equivalent 
form given by ()2.H2|1 . We set up the 1/Nc expansion for this theory as follows. For 
this theory, the graphical representation of the Feynman rules are given in Fig. |21 
Here the propagators and the vertices are easily read off from the total Lagrangian 
()2.32|) . Note that there is a cross-propagator between and Bap for p ^ 0, i.e., 
except for the another formulation. 

From the viewpoint of the l/N^ expansion, the field G^y behave as from the 
identity 

Other fields y?^, V^^ behave as S^^, in the large Nc, since the identity holds 

^ACD^BCD ^ ^AB^j^2 _ 4) /^^^ ^g ^) 

which is similar for large to the relationship 

fACDfBCD ^^AB^^_ (3.3) 

This is also the case for (ppp and (ppp and for and Therefore, the leading 
contributions in the condensed vacuum 00 7^ come from the interaction of with 
^^u^ '^^i '^^ i'^ Landau gauge A = and also from ^pp and i{>pp in the 

non-Landau gauge A 7^ 0. 

■^This Lagrangian is obtained from H2.29(l by setting p = and adding the three-gluon interaction 
term. Therefore, the modified total Lagrangian is also obtained from (|2.32() by the same procedure. 
Note that setting p = is equivalent to eliminating the cross propagator between and B^p, in 
the diagrams. 
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Figure 2: Feynman rules of the 1/Nc expansion with the 't Hooft couphng g"^. 
(a) Propagators for the gluon (wavy hne), two set of auxihary fields B^^, ip^, V 



fJ,U 



(zigzag line) and and (j), G^u (broken line), and FP ghost, antighost ^ , ^ , the aux- 
iliary fields (j)FP,^FP (solid line) (b) Vertices for two gluons and one auxiliary 
fields from a set 5^^, ip^, or a set 0, G^^, and vertices for ghost-antighost and one 
gluon field or auxiliary fields (l)pp,ippp. There are no four-point interaction vertices. 



In the following, we attempt to sum the leading order diagrams in the 1/Nc ex- 
pansion. We have rewritten the original Yang-Mills theory into an equivalent theory 
with the generating functional: 

Z[J] := J V<pVB^^VG,,Vy^^VV^iV^^V^^[V<PFpV^ip] exp {-Se} , (3.4) 

where the source terms are included for the original fields and those for the auxiliary 
fields are omitted here (but they will be introduced later). 

First, for simplicity, we consider only a sector containing (p and ^^f- The gener- 
ating functional is given by 

Z[J] := J V<Pexp{-SE[<P;J]}, (3.5) 
5.[0;7]:^/^P^0(.)^ + ^trln[^-] 

+ l^J^l^'T'J^, (3.6) 

where the gluon field has been integrated out and its contribution is replaced by 
the source term J^, and we have introduced the notation := J d^x and 

jr^40](x) := -(1 - p'){d^5,, - d,d,) - \~^d,d, + 5^,ct>. (3.7) 
Here we have made the replacement 

^ gVNc (3.8) 

with the 't Hooft coupling g'^ being fixed, and we have redefined the auxiliary field 
— s> 0/(T(^, which corresponds to the identification: 

m - -a^^<(xX(x) = -3(a2 - l)^M^{x)^^{x)]. (3.9) 
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The field </> is a single component color singlet field, whereas each gluon field s^^ := 
s^^T"^ is an adjoint field with iV^ — 1 components. The normalization factor in ()3.9|) 
provides the natural normalization, since the field produces a color singlet state 
tr[=e^(x)^(x)] with unit amplitude. Note that the first two terms in S'£;[0; J] are 
proportional to A*"^ for large Nc- The effective action S'e[0; J] has no hidden factors 
of Nc in this sector. Thus the overall factor of Nc does imply that the large Nc limit 
is the same as the semiclassical limit of an effective theory of color singlet field 
(glueball) at least in this sector. 

For the current theory to be meaningful in the limit Nc — > oo, therefore, the 
expansion around iVc = oo in the integral of ()3.5p can be performed by means of the 
steepest descent method (or saddle point method). The stationary point (poi^x) of 
is obtained from the saddle-point equation: 

3(^2 i 1)^2 ^^0(0^) + ^tr(jr^;^[0o(x)]) = 0. (3.10) 

This is nothing but the gap equation. The regularized form is given in the next section 
under the dimensional regularization. This equation has two branches, the trivial 
branch = and the nontrivial branch 7^ for arbitrary values of parameters. 
The nontrivial solution of ()3.10p agrees with the stationary point 0o away from the 
origin of the effective potential. It turns out that the nontrivial solution is preferred, 
since it lowers the energy of the system, as shown in the next section. If the vacuum 
is translation invariant, the stationary value 0o(a;) does not depend on x. 

Note that {s^^f is not a gauge invariant quantity. Nevertheless, it is possible to 
extract the gauge-invariant part from (.2^)^ which is written as the sum of a gauge- 
invariant part and a non-invariant part, which implies the decomposition 0(x) = 
4'inv{,x) + 0non(a:) is possible [Hj where 0i„,i, is the gauge invariant part and 0non is 
the gauge non- invariant, but BRST exact part. Then the vacuum expectation value 
(VEV) of the gauge non- invariant part vanishes {4>non{x)) = 0, and consequently only 
the VEV of the gauge invariant part gives a non- vanishing value {(f>{x)) = {(f>inv{x)) 7^ 
0, which is identified with 0o [T7] . 

Substituting 0(x) = 0o + ]^0(a^) into the action Se, and expanding the action 

around the stationary point 0o in powers of 0(a;), we obtain (up to field-independent 
constants) 

Se=N.I | 3(^,i^^., <^o + ^tr(jr,;M0o])}0(a:) 

+ I ^J^l^nU-'J^ + 0{1/Nc). (3.11) 

The first term of ()3.1H) shows that the saddle-point equation ()3.10p is equivalent to 
requiring cancellation of the tadpole term. See Fig. El for the graphical representation. 

The connected n-point correlation function is obtained by differentiating n times 
the generating function W[J] := lnZ'[J] with respect to the source. For example, the 
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X + 




Figure 3: Elimination of the tadpole contribution. The broken line denotes the 
auxiliary field and the wavy line the gluon field 



gluon propagator is obtained from 

In the expansion, the gluon propagator ^^yik'^^ in momentum space is read off 
from the source term in ()3.(ij) as 



^0 



^ {t. rPLik) + — , (3.13) 

where we have introduced the rescaled gauge parameter A := (1 — p^)A, the transverse 
projection -P^(p) := —VnVu/v^ and the longitudinal one P^^{p) := Pfj.Pu/p^- This 
result shows that the gluon acquires the dynamical mass given by 

M2 := (3.14) 
1 — p'^ 

It is important to remark that in the first term in the right-hand side (RHS) 
of ()3.1H1 or (j3.1()|l . there exists the logarithmic divergence coming from the second 
term |tr(,j^~^[0o]), which can be absorbed into the renormalization of the first term 
3(a-2-i)g2 • The same divergence appears in the second term |tr(.j^~^[0o]=>^^^[0o]) 
in the RHS of 1)3.111) . Therefore, a should be regarded as the renormalized one 
in what follows. See Appendix A for explicit verification of this fact based on the 
dimensional regularization in the modified minimal subtraction scheme (MS).^ 

■*Tlie quadratic divergence A^^ appearing in the term itr[..y^^^[0o]](/' with the ultraviolet mo- 
mentum cutoff A can be removed by introducing the countcrtcrm of the form. Ccf) (the same form as 
the source term J^(f) for 0). In the dimensional regularization, such quadratic divergence is hidden 
and does not appear. For this purpose, we can introduce the auxiliary field (j> in (|2.12a(l by replacing 
(j) with + C for the translation invariant measure Vcj). Then we have additional terms. ^Cx/^ + C(j) 
which is equal to ^Cs^^ + Ccj) apart form an irrelevant constant Cct>o- The redundant term ^Cjz/^ 
here is removed by introducing the counterterm of Curci-Ferrari type, —■^C^/^f, (in the Landau 

gauge A = 0) in the bare Lagrangian. In other words, removing the quadratic divergence A'^cj) is 
equivalent to introducing the counterterm ^A^^^. It should be remarked that the non-zero gluon 
mass is generated from the term ^<j>oS^^ by further removing the logarithmic divergence, irrespective 
of the quadratic divergence. 
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The connected correlation function {(f){x)] (f){y))conn shows the factorization 

{<P{x);cP{y)),onn := - = -^{4>{x)4>{y))- (3-15) 

c 

The auxihary field propagator ^^^{k"^) is read off from the second integrand of ()3.11|) 
based on the formula ()3.12|) after introducing the source term 7^0: 

%(^') = { 3(^2 ! 1)^2 -^/(0i^.^Mp;'^o]^^;Mp + fe;'^o]} . (3.16) 

Thus the auxiliary field propagator ^^^(A;^) has the well-defined UV behavior due 
to the cancellation of the UV divergence between two terms. It turns out that the 
propagator S'^^{k'^) corresponds to the summation of the infinite geometric series of 
the massive gluon-loop chain, see Fig. HI 




Figure 4: The propagator (j3.16p of the auxiliary field as the summation of the 
infinite geometric series of the massive gluon-loop chain. 



In the small momentum k"^ region, the propagator for A = or 1 is calculated in 
the form: 

_ 192.^(3 + A)-H1-pTM^ 

by using the gap equation ()3.1()j) and the result of Appendix A. The connected cor- 
relation function 0(?/))conn decays exponentially for large \x — y\, since the 
right-hand-side of ()3.15p is equal to the Fourier transform of ()3.17|) . Here it should 
be understood that the contributions from the gauge non-invariant part vanish in the 
final expression. In this way, we can obtain the mass of the lightest scalar glueball 
O"^"*". The ratio of the glueball mass 0"^"'" to the gluon mass is given by 

M(0++)/M = ^ 2.45, (3.18) 



in the leading order. This ratio is gauge independent, i.e., independent of the gauge 
fixing parameter A and it does not depend on the parameter a. This value is similar 
to the prediction of the potential model of Cornwall and Soni |39|l4nj : M{0~^~^)/M = 
2.3. If one adopts the gluon mass M = 600 ~ TOOMeV as suggested from the 
numerical simulation on a lattice [H] , this ratio reproduces the expected value of the 
glueball mass jl2] M{0~^~^) = 1.5 ~ 1.7 GeV. Heavier gluon mass around 1 GeV 
is also suggested, see phenomenological analyses (121; lattice simulation [33] and the 
analytical study 
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Next, we include the antisymmetric tensor field B^^^ and consider a sector consist- 



ing of 0, B^j^ and s^/^. The generating functional is given by 

Z[J] ■.= IVcl)VB^,exp{-SE}, 
Se := I -^B^,{x)-B^^%x) + 



6(^2-1)^2 



AB-i 



(3.19) 



where 



AB 



9 r 



ABC*j^C 



(3.20) 



The action has the 1/Nc expansion around 0o (up to field-independent constants): 



S'e=N, 



3(a2 - 1)^ 



tr(jr^;i[0o]) 



3((t2 



1 



'B 



*B 



utB 



(3.21) 



NcS"^^. In (KT^ . a contribution coming from 



where we have defined I/j^uais '■= hi^paSuis — ^pp^ua), and we have used the identities: 
^AB^AB = Nl-l and /^'^^/^ " 
the term is simplified as 

1 



-p-d>:B^^^^.J^-^[M5^^dr*B 



utB 



2-(l-p2)92 + 0^ 



(^9^*^^""^, (3.22) 



where the gauge-parameter dependent part in ^^^[(f>o] disappears due to the anti- 
symmetric property of B^,^. Therefore, the kinetic term for B^,^ is generated and the 
auxiliary field becomes dynamical in the leading order oi l/N^ expansion irrespective 
of the magnitude of the 't Hooft coupling g, as far as 0o 7^ 0.^ 

In the action, the terms with more than three auxiliary fields are suppressed in 
l/Nc- Hence the leading contribution is exhausted by a few terms. However, such 
simplicity is apprarent and is immediately lost, as soon as we proceed to obtain the 
propagator. In fact, the auxiliary tensor field propagator Bq^(^^) the leading 

^This phenomenon has already been claimed to occur in |25II26| based on the perturbation ex- 
pansion in the coupling constant to a certain order. Therefore, it was unreliable to extend the result 
to the infrared region where the coupling constant becomes quite large. 
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order consists of all the planar diagrams with two external lines corresponding to 
and Baj3 (see Fig. El) where the bare propagator is given by^ 

and the vertex is read off from the action ()3.19|) and ()3.20|) . Such planar diagrams 
are obtained from the gluon one-loop diagram with a number of external fields B 
by contracting the external lines two by two. In each step, a factor N is increased 
without increasing the order of g. All the diagrams in Fig. are the same order in 
the 1/Nc expansion as can be seen by making use of the double line representation. 

The similar remark is also applied to the propagator of the auxiliary field (p. The 
inverse propagator of the auxiliary field </> in the leading order is the sum of the similar 
diagrams to Fig. El where two external B lines are replaced by the 0. This is because 
the graphs with internal </)-lines give higher oder contributions in l/Nc- 




+ 



Figure 5: Planar diagrams in the leading order for the inverse propagator of the 
auxiliary tensor field 5^,^ where the zig-zag line denotes the S^jy field and the wavy 
line the gluon field. 

Thus, it is still impossible to sum up all the planar diagrams in question to have 
the closed form for the propagator in the leading of the 1/Nc expansion. Here we 
obtain the auxiliary tensor field propagator, the inverse of which is truncated up to 
0{g^) as above 

^Znjk') := I^W/^^^" + p + W/^^^" + KUk)} \ (3-25) 

^We have used the integration by parts 

= I - B^^d.d.B^^ + B^^d.dxBt). (3.23) 

For simphcity, we have imposed the gauge-fixing condition d\B\u = 0, although this is not an 
indispensable condition, see e.g., |26| for more details. 
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where Yl^^^^^k) is the vacuum polarization tensor in the lowest order given by 



•'0j-'a/3,p2O'2 



(3.26) 



The explicit form of the vacuum polarization tensor n^jf„^(/c) is given based on MS 
in Appendix B. The obtained propagator corresponds to the summation of the in- 
finite geometric series of the massive gluon-loop chain, see Fig. IHl for the graphical 
representation. This propagator obtained by the partial resummation of the lead- 
ing order graphs is called the truncated propagator hearafter. It is possible to resum 
the diagrams based on the Schwinger- Dyson equation approach which will be given 
elsewhere 




_|_ A/Vg^^|^~^^ 



Figure 6: The propagator ()3.25|) of the auxiliary tensor field B^^i, as the summation 
of the infinite geometric series of the massive gluon-loop chain. Here the sum of the 
first two propagators denotes the propagator ()3.24|1 . 



In the small momentum region, the vacuum polarization tensor for Bf^y reads^ 



{k) = -6 



AB 9 



a 



167r2 (1 -p2)2 

M6 J ' 



<■ fiu,af3 



1 + Al P k^ 



(3.27) 



where /o = ^ In 



Af2 



4' 



/l 



11 

180' 



/2 



105 



for the Landau gauge A = and 



/o = In /i = /2 = yHo Feynman gauge A = 1. Thus, the truncated 

propagator ^b^,b^^(^^) of the auxiliary field has the IR asymptotic behavior 

r 1 ^2 rr2 1"-^ 



2((t2-1) 167r2(l-p2)2 



fl 



M2 



m2 + P + 0{k^ 



where 



m 



1 



9' 



a 



2(a2 - 1) 16712 (1 - p2)2 



/i 



1 



9' 



2 16712 (l-p2)2 



/o 



(3.28) 



(3.29) 



^The divergent term is absorbed into the rcnormalization of the tensor field B^^ through the 
renormahzation constant Zb defined by B^i, = Z^J^B^^, just as in the case of Maximal Abelian 

gauge mmni. 
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In the large momentum k"^ region, on the other hand, the truncated propagator 
reads 



(3.30) 

If such a high energy behavior as ()3.3()p is preserved after the inclusion of all the 
planar diagrams, our theory of gluons with non-vanishing effective mass is expected 
to be renormalizable order by order in the 1/Nc expansion. This is in sharp contrast 
to the naive massive Yang-Mills theory which is known to be nonrenormalizable in 
weak-coupling perturbation theory. More arguments on this point will be given in a 
subsequent paper. 

Other fields can be treated in the similar way. In fact, after an appropriate 
rescaling of other auxiliary fields as (j), the total action reads 

^EYM = I [\b^, ■ B^'^ + ^ IS^G^^^G^- 



+ ^ ■ 9^(9^^ + y/a^(j)Fp'>^ ■ ^ + ^(t'fp^fp ■ (^ * ^) 

1 ^Ar 1 



2 



/M^.'']-'/.''} + 2 l^Det[jr,f^], (3.3i; 



where 



^AB . ^ \A\X.AB ■9_„eABC*nC , AT-^xABx 1 , /T^MBC j: ,„C 



■■=Jt - wdu^B,''^ + (1 - i)^{d,^ X ^)^ - X d,^)^. (3.32) 

4 Estimate of gluon mass and vacuum condensate 

Now, we examine based on the effective potential whether the vacuum condensate 
of mass dimension two is realized or not. All the auxiliary fields B^^, (p, ip^, G^y, 
have the canonical mass dimension 2. Therefore, various vacuum condensates 
can contribute to lower the energy of the system. However, if the color rotation 
invariance and Lorentz rotation invariance are required, we find {B^^) = {(p^) = 
{Gnu) = {V^) = 0. Only a possibility examined in this paper is (0) 0. This 
vacuum condensate corresponds to {{^ff'Y) 7^ which breaks neither color symmetry 
nor Lorentz symmetry. In what follows, we treat the vacuum condensate of mass 
dimension 2 alone.* 



^This does not exclude the dimension four condensate of the type: (B^^B^^), {ip ip ), (G^^G^iy), 
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The minimum value of the effective potential should not depend on the gauge 
fixing parameter. Therefore, we calculate only the A = case. In this case, the 
BRST-invariant composite operator reduces to .e^. 

The effective potential is obtained from^ 



V 



2 ^ 



■lndet[-(l-p^)5^^(9^5^,-9^9, 



■ (4.1) 



The scalar part of the effective potential is obtained by resumming the gluon one-loop 
contribution as 



2 ^ 2 



(2^ 
(2^ 



In det 



(1 - p'){pH>''' - p^pn + jp'^p" + 



In 



(1 - p^)p' + p' + Xa^(p /X 



(4.2) 



At A = (up to (^-independent constant), we obtain 



=^^2 , 3(iV^^ - 1) 
where V{0) = 0. The first derivative reads 



d'^p 

1 -p2 



In 



In 



l-p2 



47r/i2 



(4.3) 



V'{<t>) =0 



3(iV< 



647r2 



327r2 



In 



47r/i2 



(4.4) 



The stationary condition V^'((/>) = has two solutions, = or = 
Fig. 13 Thus the absolute minimum of the potential is given at — 0o 



/>o 7^ 0. See 



1 — p^) =47r/i^ exp 



=47r/i exp 



327r2 (1 - p2)2 



3(iV2-l)a^ 2 J 

iV2 (0-2 - 1) 327r2 _ 1 
9o2 ~ 2 



(iV,2 



(4.5) 



Note that the right-hand-side has a finite value in the limit N^. oo. It turns out 
that 0"^ must be larger than one, o"^ > 1. This result can be compared with the 
perturbative calculation of the effective potential to the three- loop order [IH] . 

Now we estimate the order of the gluon mass by combining the above results with 
the lattice data. Using the beta function at one-loop level in the relationship of the 
trace anomaly 



4a., 



4q;? 



(0|T/|0) =^mKf\0) = ^^(0|f (^,1.)^|0) 



— ^(o|^(^;L)», 



(4.6) 



^We use the same letters to denote the field obtained after the Legendre transformation. 
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and the values of the gluon condensation of mass dimension 4 calculated on a lat- 
tice mi 

(0|-(=^„1.)'|0) = 0.152(GeV)'(iV, = 2), 0.144(GeV)'(iV, = 3) (4.7) 

we obtain 

(0|T/|0) = -0.1392(GeV)^(A/'c = 2), -0.1980(GeV)^(iVe = 3) (4.8) 

where we have put the flavor number A'^^ = 0. 

The effective potential just obtained has the minimum at — 0o- 

1287r2 \l-p^J 1287r2 ^ ^ ^ 

Supposing all the contribution to the vacuum energy comes only from the condensate 
00, i-e., 



V{<Po) = ^(0|T/|0) = -0.0495(GeV)'(for = 3), (4.10) 



we obtain the rough estimate of the effective gluon mass 



= V^</,0o/(l - ^-2) = 1-27 GeV. (4.11) 

This value does not depend on the specific choice of the parameters p and a. The 
estimated value of the effective gluon mass is larger than the result: rrij^ = 0.6 GeV 
of numerical simulation on a lattice [H] . The vacuum condensate of mass dimension 
2 is given by 

K/OI =2/00 = = ^m^. (4.12) 

0-0 2cr^ 
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5 Static potential 

We introduce the source J^^ of B^y by modifying the identity ()2.4|1 into 

1 = j VB^, exp - j d^x^{B^, - i*[p{d^K - dus^^) + ag{£^^ x <)] + J^.f 

(5.1) 

This leads to the source term 

\b,,J^- - '-*r-p{d,K - dys^,,) + \{J,u)\ (5.2) 



and the modification of Jf^^ into Jf^^^ X^J^ - tagf^^^*J^''^. 

In the similar way, we can introduce the source J"^, J'^f_iu, J"^"^, J^^^ auxil- 



iary field 0, G^^, ip^, V^. Then the source term is further modified by adding 



+ ■ + -X. ■ + T^iJ]^^)"- (5-3) 



Thus, the introduction of the source terms for the auxiliary fields is equivalent to the 
replacement 

and the addition of the terms 



The linear term in the gluon field appears only in the source term ()5.2j) of B^^, 
which is recast into the form 

J d^x^B^, ■J'^^ - j ^"^-^ ■ J'^^ = -^P^''*J^.v. (5.6) 

which cancels the original source term / d^x-c^ ■ for s^^. The current defined 
by ()5.6|1 satisfies the conservation law 

d^J^ = 0. (5.7) 
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due to antisymmetry of the tensor current, J^^, 
()5.7|) . the J^i, can be solved in terms of J^: 



~Ji/ti- With the conservation law 



J, 



where ra^ is a constant vector corresponding to a set of integration constants. There- 
fore, we can calculate the generating functional for the source for the original gluon 
field ^ by making instead use of the antisymmetric tensor current source J^^ cou- 
pled to the antisymmetric tensor field B^y via the relationship between two currents 

(lEHD. 

The integration over the field B^y is performed as follows. 

/ 1)5,. exp {- / d'x [^*B^'^\&^l^^J-'*B^^^ + \*B^^:r'^^ 



exp 
exp 

exp 



where r 



1 2?7^ 



d^k -1 

(27r)4 4p2 



Jf. 



^B,B{k)- 



n 



n ■ k) 



2 M 



(5.9) 




<l-> 



t 
A 




>- X 



(a) 



(b) 



Figure 8: Two equivalent descriptions for calculating the static potential, (a) J'^ 
coupled to £/fj,, (b) J^'^ coupled to Bf^y. 



We consider a pair of opposite color charges and —Q^ located at = a 
and X = b, respectively. In this case, the conserved current is given by J'^i^x) = 
Q'^6^q[6'^{x — a) — 6^{x — b)]. Using the Fourier transform 

J^{k) = Q^5,o27r5(A;°)(e-'^'^ - e^^'^"), (5.10) 

we obtain 



1 r r d^k 1 
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where r := a — b is the relative vector between a pair of charges, and we have chosen 
n'^ = (0, n) with a unit vector n: \n\ = 1. Therefore the static potential is calculated 
from the formula: 



(0,fc)- 



n ■ ky 



:i-e 



-ik-r\ 



where we have used ^s,s(0, —k) = ^^b,b{0, k). 

Substituting for ^b,b{0, k) the = hmit of the IR propagator (jS^Hl): 



(5.12) 



2(a2-l) 167r2 (1 - f7-2)2 



fi 



M^{m' + e + 0{k^)} ' 



a 



n -1 r 



2(^2-1) 167r2 (1 - (7-2)2 

:=0o/(l-^-'), 



A 



1 



2 167r2 (1 - (7-2) 



(5.13) 



into the above equation ()5.12j) . we obtain the linear static potential as a function of 
r := \r\ for r ^ 1^° 



Vdr) 



(5.14) 



where the string tension (Jst (with an enhance factor) is given by 



(^st =0" 



-a 



fcr2 - I] 



1 - 



87r2 (a2 



^7r2 (1 
-/i 



1 




(5.15) 



with A being the ultraviolet cutoff. In the classical limit a | 1, the string ten- 
sion vanishes agt [ 0. The string tension is essentially proportional to the vacuum 
condensate 0o. Thus, the non-vanishing string tension is generated from quantum 
effects. 

This string tension has the very similar form to that obtained by making use of 
Abelian projection assuming the Abelian dominance [^EHj up to a numerical factor, 
as expected from the type-II dual superconductor. ^2 xhe constant vector was 

^"The vector n is chosen to be parallel to r by the axial symmetry of the system in question and 
the minimum energy condition. The method of calculations in what follows is the same as those 
of |48[I49| . See there references for the details of calculations. 

-'^^This UV divergence originates from the fact that the VEV 4iq{x) = {4>{x)) is a x-independent 
constant. In other words, the world sheet of the effective string connecting two color charges has 
no width. The presence of the UV cutoff A means that the string has a non-zero width of the order 
A~^. See the next section. 

^^It should be remarked that we have used the Abelian-like expression for the current for 
defining the Wilson loop. To incorporate the full non-Abelian structure, we must use the non- 
Abelian Stokes theorem which enables us to rewrite the line integral in the non-Abelian Wilson loop 
into the surface integral over the surface whose boundary coincides with the Wilson loop. Such a 
calculation was performed in the previous paper |2f>j . In this case, only the maximal torus subgroup 
is responsible for the final result, even if we start from the Lorentz gauge. Thus we reproduce the 
dual superconductor picture. 
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introduced following the Zwanziger formulation which allows the simultaneous 
existence of electric and magnetic charge. Under the Abelian projection, only the 
Abelian charges corresponding to the maximal torus group contribute to the string 
tension. In the Lorentz gauge, however, all the components of color charges Q"^ 
contribute to the string tension. Thus the propagation of B^y [Boi in this case) gives 
rise to the linear potential. Incidentally, in order to know the potential Viir) for 
small r, we must use the UV asymptotic propagator of the form 

On the other hand, the original source term for the gluon field ^ 

j4^~-'- jA\ \i ^<:AB^-l tB f tAi /I „2\q2 , i 1-1 tA 



yields the Yukawa- type potential (with an enhance factor due to a): 

Vy{r) = -(l-a-r'^^, (5.17) 

Air r 

apart from an irrelevant constant. The Yukawa-type potential does not depend on 
n'^ and it is gauge-parameter independent due to the conservation law ()5.7|) . 

Here we comment the electric-magnetic dual picture of the Yang-Mills vacuum. 
The dual transformation (or Fourier transformation) can be used to transform the 
massive antisymmetric tensor field theory into another equivalent theory. It is well 
known that in -D-dimensional space-time, a massless antisymmetric tensor field of 
rankp (p-form) is dual to the [D — p — 2)-foTm, while a massive antisymmetric tensor 
field of rank p is dual to the [D — p — l)-form. In fact, it is shown that the effective 
theory retained up to the quadratic terms in the massive Kalb-Ramond [^21 field B^^ 
is dual to the massive vector field f/^ theory described by the dual Ginzburg-Landau 
(DGL) theory in the London limit. See [23123123 for details. Therefore, the quite 
similar situation to the dual superconductivity is realized in the Yang-Mills vacuum in 
the 1/Nc expansion and the the London limit corresponds to in the leading order of the 
expansion. To go beyond the London limit, it it necessary |25tl26] to incorporate the 
quartic term in B^^ which appears in the next-to-leading order of the 1/iVc expansion. 



6 String representation 

In the derivation of static potential in the previous section, we adopted a specific 
configuration for the color charges and introduced a constant vector n^. These short- 
coming are avoided and we can derive more general result by means of the string 
representation of Yang- Mills theory as will be shown in what follows. 

The trajectory C of the charge is parameterized by a parameter r. Hence the 
current is written as 

J^ix) = Q^l^dy,iT)S\x-yir)) = j^dr^^5\x - y{r)). (6.1) 

Let the trajectory be a loop C and let S be the surface whose boundary is C, i.e., 
9S = C . The surface S is parameterized by two parameters r, a. Then it is shown 
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that the vorticity tensor current defined by 



drda 



satisfies the relationship 



(6.2) 



(6.3) 



The Q^y{x) has its support on the two-dimensional surface E, while the J^(x) has 
the support on the loop C in the four-dimensional target space x^. As in the previous 
section, we have the correspondence 

(6.4) 



The contribution of the vorticity current is estimated as 
/ Pi?,.exp {- j d'x [-*B^^\&^l^^J-^*B-^^ H 
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X 
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(6.5) 



The action in ()6.5|1 describes the interaction between two surface-elements dS^y 
and dSaii in the surface S spanned by the loop C . This expression fl6.5|l is the 
same as the non-local action for the vorticity tensor current, from which the string 
representation of the Wilson loop average in Yang-Mills theory was derived in 
by using the method j2ZI- We can show 

(jnUl) =exp(-S'cs[x]), 



Scs =(^st d'^a^/g + aQ^ d^a^/gg^'^dat^.^dbt^^ ^ '^^ X ^^^^^ ' ^^'^^"^ 



where 
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fi 



/>o In 1 + — 



-g^(i-a-Y<o, 



(6.6b) 

(6.6c) 
(6.6d) 



In the action ()6.6a|l . the first term is the Nambu-Goto action, the second term is 
the rigidity term, and the third term is the intrinsic curvature term. See 
for the details of the calculation. The string theory just derived is nothing but the 
rigid string with a negative rigidity term proposed by Polyakov [2H1- Thus we have 
shown that the vacuum condensate of mass dimension two gives a contribution to the 
non-vanishing string tension (Jst in the linear potential already in the leading order of 
expansion. 
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7 Conclusion and discussion 



We have proposed an equivalent formulation of the Yang-Mills theory in which all 
the gluon self-interactions are absent. In other words, the gluon fields are at most 
quadratic in the new formulation and all the gluon fields can be integrated out exactly 
in the SU{Nc) Yang-Mills theory to obtain an equivalent formulation at the cost of 
introducing various auxiliary fields, i.e., an antisymmetric tensor -B^j,, two symmetric 
tensors G^u, V^, four Lorentz scalars 0, ip^, 4>fp, '^fp- Hence the resulting theory is 
written in terms of the seven auxiliary fields and the ghost ^ and antighost fields 
Two of the auxihary fields 4>FP,fFP necessary only when the four- ghost 
interaction exists. 

Then we have set up a new l/N^. expansion on a non-trivial vacuum in which 
we have shown that one of the scalar field (p can acquire a non-vanishing vacuum 
expectation value {(p) = (pQ ^ 0. This corresponds to the on-shell BRST-invariant 
vacuum condensate of mass dimension two pHlElE] which reduces in the Landau 
gauge A = to the gluon pair condensation 7^ 0. Once such type of Bose- 

Einstein condensation occurs, all the gluons acquire a common effective mass M. 
This behavior of the gluon propagator is consistent with the recent results of the 
Schwinger-Dyson equation |^ and the axiomatic approach |55j . 

In the condensed vacuum, all the auxiliary fields acquire their kinetic term and 
consequently become dynamical fields. We have given the Feynman rules of the l/N^. 
expansion which give the propagators of gluon fields, auxiliary fields and (anti)ghost 
fields, and the interaction vertices among the gluon and other fields. The resulting 
theory is expected to be renormalizable in the scheme of the 1/Nc expansion just 
obtained. However, the detailed analyses of the renormalizability will be postponed 
to a subsequent paper. 

The effective gluon mass has been estimated in the leading order by assuming that 
the vacuum energy is saturated by the vacuum condensate (f>o alone. The obtained 
value is larger than the expected value due to numerical simulations on a lattice. This 
might suggest that we must go beyond the leading order or take into account other 
vacuum condensates. However, the result does not contradict with the experimental 
upper bound Moreover, the renormalization of the composite operator in 
the given 1/Nc expansion is to be investigated in future works. 

As a by-product, the mass of the lightest scalar glueball mass to the gluon mass 
is obtained by performing the partial summation of the leading order diagrams of the 
1/Nc expansion: M{0~^^)/M = y/E = 2.45 independently of the parameter a. This 
value agrees with that of the potential model by Cornwall and Soni [2211101 • In order 
to reproduce the expected mass of the glueball 0"*"+, i.e., M{0^^) = 1.5 ~ 1.7GeV, the 
gluon must acquire the mass M = 600 ~ 700MeV in agreement with the numerical 
simulation on a lattice [H]. It is desirable to calculate the complete glueball spectrum 
in the future work. 

The static quark potential calculated in the same approximation is expressed 
as a combination of the Yukawa type potential and the linear potential. The non- 
vanishing string tension is provided by the vacuum condensate of mass dimension 
two proposed by ^U]- It is argued that the mixed gluon-ghost condensate as the 
VEV of the on-shell BRST invariant composite operator of mass dimension two JU] 
means the transverse gluon pair condensation of mass dimension 2, which is gauge- 
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invariant but nonlocal ^2)- In the Maximal Abelian gauge, these results are consistent 
with the dual superconductor picture of QCD vacuum examined through the string 
representation 

The obtained string tension ast depends on the parameter a, just as the gluon 
mass M does. This result should be compared with the other large results, e.g., 
the N = oo extrapolation of the SU(2) and SU(3) simulation results on the lat- 
tice ISE]: M{0^^)/y^a^ = 3.56 ± 0.18, ■ ■ ■ , or the light-front Hamiltonian approach 
of the transverse lattice gauge theory [SZj: M{0~^~^)/ ^/a^ = 3.3 ± 0.2, ■ ■ ■ . In the 
phenomenological applications, the parameter a and A (in the string tension) can be 
adjusted to reproduce the gluon mass M = 600 ~ 700MeV and the string tension 
^/a^ = 450MeV. By adopting these values as inputs, one can calculate the complete 
glueball spectrum and other physical quantities. 

It is also interesting to study the fate of the pair condensation in relation with 
deconfinement , if the sufficient number Nf of quark flavors are incorporated into the 
theory. It is interesting to confirm whether or not the pair condensation disappears 
above a critical number of quark fiavors jS]. Some phenomenological consequences 
due to the vacuum condensate of mass dimension two were discussed also in the 
papers [SHj . 

Note added: After having submitted this paper to e-Print archive, the author 
was informed that some similar results were previously obtained on the relation among 
the condensate of dimension 2, the gluon condensate and the mass of the dual gluon, 
which is related to the glueball jll] (but just inside two models, i.e., stochastic vac- 
uum model inni and dual long-distance QCD [00], and not with a complete theory 
derivation as this paper). The author would like to thank Nora Brambilla for sending 
this information. 



A Cancellation of divergence and renormalization 



First, we focus on the divergence in 

1 



3(^2 _ 1)^2- ■ 2 

The dimensional regularization yields 

d^k 1 r(a-f) 



2 0o + ^tr(jr^;i[0o]). 



«+D _ r(e-2 + a)^^ ^2_„_^ 



(27r)^ + $0]" (47r)^/2r(a)'^"' (47r)2-"-^ "^"^^^ 

where e := 2 — D/2. Defining $0 '■= 0o/(l — P^) and A := (1 — p^)A, we have 
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where we have defined e ^:=e ^ — 7 + ln(47r) and used the Laurent expansion of the 
Gamma function r(e — 1) = — e^^ + 7 — 1 + 0(e). Then we obtain 



1 



+ (D - 1 + P) In 



D - 1 + 



l-p2 



1 +PlnA + 0(e) 



(A.4) 



Therefore, the logarithmic divergence of the second term should be absorbed into the 
first term, 



divergent part of 



1 
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Next, in order to calculate the divergent part in the second term of 



1 



3(a2-l)^2 2 J (27r)^ '^^ 
we put k = 0. Then, using the integration formula ()A.2|) . the second term reads 
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where we have used r(e) = e ^ — 7 + 0(e). Therefore, the logarithmic divergence of 



the second term is cancelled with that of the first term ^ 
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provided that the relation ()A.5jl is satisfied. Therefore, the auxiliary field propagator 
^4>4>(^'^) well-defined finite UV expression. 

In order to absorb the divergent part in the relation ()A.5|) . we perform the renor- 
malization. We introduce the renormalization factor Z„ hj a = Z^aji to absorb the 
divergence in such a way that 



divergent part of 
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the renormalization factor is determined as 



,,,__3^r^ D_ 
~ 4 (47r)2 ^« 



(A.ll) 



where we have introduced the renormahzation scale [i to render the couphng constant 
dimensionless for arbitrary D. The RG function is calculated as 
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This implies that increases as n decreases. This tendency is favorable to our choice 
cr^ > 1. The explicit solution ctr of this equation for D = 4 and A = is given as 



(Tij(/i) = 0-i?(/io) 



1 + a-ij(/io 
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(A.13) 



by making use of the one- loop expression of the running coupling constant JTj. This 
equation is valid when A ^ < /iq, i.e., for small 't Hooft coupling g"^. Note that 
the Nf. dependence enters in this expression through the 't Hooft coupling g"^. 

At a certain value of /i, say, /xq, we can choose an arbitrary value of a^, say, 
c"o '■= d'Rif^o)- Then the value of c3"^(/i) at another point fi is calculated according to 
the RG equation. Here the initial value Pr{hq) of Pr^h) at /i = /iq must be tuned so 
that Pr{h) = cr]^{jji) at the desired value of since the gluon self-interactions are 
eliminated only when this relation is satisfied. In the 1/Nc expansion, the magnitude 
of 't Hooft coupling g^ can not be specified, although it is expected to be small in our 
treatment. 



B Calculation of loop integrals 

The dimensional regularization in the modified minimal subtraction scheme (MS) 
(e~^ := 2/(4 — D) — '^e + ln47r) leads to the following results in the Euclidean for- 
mulation. In following equations, fc^ denotes the Euclidean momentum squared. The 
corresponding expression in the Minkowski region is obtained by k'^ —>■ —k"^ below the 
threshold k"^ < 4M^. It is possible to obtain the full expression for arbitrary A. But it 
is quite complicated and long. Hence we write down two special but important cases 
explicitly. 

(1) In the Landau gauge A = 0: 
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where := 0o/(l — p^) and we have introduced 

p }l(pT pT _ pT pT \ pT e _ ^/J-^i^ 

In the large momentum k"^ region, on the other hand, the propagator reads 
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Here we have used the following properties of the projection operator I^u,ai3- 

k'^^ Ifj,u,piijk^^ Iaf3,p2<J —'^{k^ka^iip k^kp5ya kyka^^p + ki^kf^6fj^a) 
Ipu^apB'^'^ =Bfj_^, B^'^I^jj^aP = -Bo/3, 

for an arbitrary antisymmetric tensor 5 
(2) In the Feynman gauge A = 1: 
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In the large momentum k"^ region, on the other hand, the propagator reads 
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For more details of the calculations, see 
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